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ABSTRACT 


In many polarimetry applications, including observations in the X-ray band, 
the measurement of a polarization signal can be reduced to the detection and 
quantification of a deviation from uniformity of a distribution of measured an- 
gles of the form A + B cos 2 (0 — 0 O ) (0 < 0 < tt). We explore the statistics 
of such polarization measurements using both Monte Carlo simulations as well 
as analytic calculations based on the appropriate probability distributions. We 
derive relations for the number of counts required to reach a given detection 
level (parameterized by f3 the “number of a's v of the measurement) appropriate 
for measuring the modulation amplitude a by itself (single interesting parameter 
case) or jointly with the position angle 0 (two interesting parameters case). We 
show that for the former case when the intrinsic amplitude is equal to the well 
known minimum detectable polarization (MDP) it is, on average, detected at the 
3 a level. For the latter case, when one requires a joint measurement at the same 
confidence level, then more counts are needed, by a factor of ~ 2.2, than that 
required to achieve the MDP level. We find that the position angle uncertainty 
at la confidence is well described by the relation a ^ = 28.5(degrees)//3. 


Subject headings: polarimetry: general — statistical analysis: Monte Carlo simulations 



1. Introduction 


Emission and scattering processes thought to be important in many astrophysical 
X-ray sources are likely to impart specific polarization signatures, but to date there have 
only been a few positive detections of polarization from cosmic X-ray sources, largely due 
to sensitivity limitations. Some of the earliest and highest precision measurements were 
made with the OSO-8 Bragg reflection polarimeter (Kestenbaum et ah 1976; Weisskopf et 
al. 1976), and include high significance measurements of the linear polarization properties 
of the Crab nebula in several energy bands (Weisskopf et ah 1978). 

More recently, observations with the INTEGRAL spectrometer SPI and imager IBIS 
have exploited the polarization dependence of Compton scattering to infer the linear 
polarization properties of the Crab at y-ray energies (Dean et al. 2008; Forot et al. 2009). 
These results indicate that the > 200 keV flux from the Crab nebula is highly polarized 
(~ 50%), with a position angle consistent with the pulsar rotation axis. 

In the last few years the development of micropattern gas detectors has enabled the 
capability to directly image the charge tracks produced by photoelectrons, thus enabling use 
of photoelectric absorption in a detection gas as a direct probe of X-ray polarization (Costa 
et al. 2001; Black et al. 2004; 2010). It is likely that such technology will be employed in 
the not-too-distant future to sensitively explore the polarization properties of many classes 
of astrophysical X-ray sources for the Erst time. 

In this paper we explore the question of how one detects, measures, and characterizes 
a polarization signal with a photoelectric polarimeter. The remainder of this paper is 
organized as follows; in §2 we outline the basic problem of detecting a modulation in 
a distribution of angles, and we describe the angle distributions used throughout the 
paper. In §3 we briefly outline the probability distribution relevant to such polarization 
measurements. In §4 we describe our Monte Carlo simulations and present our results. We 



conclude with a short summary in §5. 


2. Statement of the Problem 

The angular distribution of photoelectrons ejected by a linearly polarized beam of 
photons (X-rays) is proportional to sin 2 (d) cos 2 (0)/(l — 0 cos(0)) 4 (see, for example, Costa et 
al. 2001 ), where 9 is the emission angle measured from the direction of the incident photon 
(0 < 9 < 7r), and 0 is the azimuthal angle measured relative to the polarization vector 
of the incident photon (see Figure 1 for the basic geometry applicable to a photoelectric 
polarimeter). In most practical situations the angle 9 is not inferred directly, but the 
electron charge track projected into the plane orthogonal to the direction of the photon 
(the plane defined by 9 = 90 deg) is imaged and thus 0 can be estimated for each detected 
photon. The angle 0 is measured around the line of sight to the target of interest and 
can be referenced to, for example, local North on the sky. The presence of a significant 
linear polarization component is then evident in the distribution of azimuthal angles. For 
example, an unpolarized photon flux will produce a uniform distribution in the angle 0, 
whereas a linear polarization component produces a distribution peaked at a particular 
azimuthal angle, 0o- 

In general, the observed distribution in 0 will be of the form, 

S(<f>) — A + B cos 2 (0 - 0 O ) . (1) 

This distribution has an amplitude of modulation given by, a = (S max — S m i n ) / (S max + 
S mtn ) = B/(2A + B ), and a position angle (the angle at which the distribution has a 
maximum) given by 0o, thus, the detection of polarization can be reduced to a statistical 
detection of a modulation in the distribution of angles, 0. Such a distribution is often 
referred to as a modulation curve. 
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In principle the angle 0 can be measured in the range from 0 — 2n (0 - 360 degrees), 
however, due to the two-fold symmetry of the cos 2 (0) dependence of the angular distribution 
of ejected photoelectrons, it is sufficient to consider distributions over a range of angles 
from 0 to 7T (0 - 180 degrees). 

An equivalent and often convenient way to express this distribution is using the 
so-called Stokes decomposition, 

S((j)) = / + Q sin {2(f)) + U cos (20) , (2) 

where /, Q and U are the well-known Stokes parameters, the modulation is now given by 
a = ( Q 2 + f/ 2 ) 1 / 2 //, and the position angle is 0 O = 1/2 taa~ 1 {Q /U). 

Now, the amplitude of modulation a is not in general equal to the source polarization 
amplitude, a p , because a detector is not a perfect analyzer and will not provide an exact 
measurement of the true photoelectron angle 0. That is, individual angle estimates will have 
some uncertainty associated with them and these will produce a uniform (unmodulated) 
component to the measured distribution even in the case of a 100% polarized beam. This 
“lossiness” of the angle estimates is quantified in terms of the so-called detector modulation, 
H , which is the modulation amplitude produced in the detector by a 100% polarized X-ray 
beam. 

In general, // can depend on a number of factors, including the energy of the incident 
photons, and the composition and pressure of the absorbing gas, among others (see 
Paciani et al. 2003, Bellazini et al. 2003). In the absence of background, the amplitude 
of polarization is just a p = a/fx. In general, a p is larger than the measured amplitude of 
modulation, a, because as noted above, detectors are not 100% efficient, and some of the 
intrinsic polarization amplitude is smeared out. 
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3. Probability Distribution 

Previous studies have shown that the joint probability distribution for a measurement 
of polarization amplitude, a, and position angle, 0 is given by 

P(N,a,a 0 ,<j>,<t>o ) = ^ exp ^-^(a 2 + «o “ 2aa o cos2((0 - 0 O )))^ , (3) 

where a, ao, 0, <j>o, and N are the measured amplitude, the true amplitude, the measured 
position angle, the true position angle, and the number of detected photons, respectively 
(see, for example, Weisskopf, Eisner & O’Dell 2010). In the case of no intrinsic polarization, 
ao = 0, the distribution simplifies substantially to 

P ( N ’ a ) = ex P > ( 4 ) 

and this expression can be readily integrated to find the probability of measuring an 
amplitude, a, if there is no intrinsic polarization. The amplitude that has a 1% chance 
of being measured is referred to as the minimum detectable amplitude (MDA), and it is 
relatively straightforward to show that MDA = 4.29 /yfN. The polarization amplitude 
that would produce this modulation amplitude in a particular detector system is called 
the minimum detectable polarization (MDP), and based on the discussion above is just 
MDP = MDA/ p = 4.29 /(fiy/N). 

Furthermore, if we are not concerned with the position angle, 0, then we can integrate 
equation (3) over angles and obtain the distribution, 

P(N , a, ao) = f P(a, <j>)d</> = ^exp(- — + 9 Q °’ > )/ 0 (^) (5) 

J-n/2 & 2 o'- a 2 

where Iq is the modified Bessel function of order zero, and cr 2 = 2/N. This distribution is 
known as the Rice distribution (Rice 1945) and it reflects the fact that the amplitude is 
always a positive quantity, and so the distribution must go to zero for a = 0. An important 
property of this distribution concerns the second moment, which is related to the width, 
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and is given by: < a 2 >= + 4/N which shows that the distribution width increases with 

a o- 


4. Monte Carlo Simulations 

Without loss of generality an X-ray source’s linear polarization characteristics can be 
described by the modulation amplitude, a 0 , and position angle, </> 0 (in the range from 0 - 
7r), that would be produced in a particular detector system. In making an observation a 
total of N photons are observed for each of which an angle </>, is estimated for the ejected 
photoclectron in the range from 0 — n (0 - 180 deg). One can then create a histogram of 
the number of events (photons) in each of M position angle bins, where the bin size (in 
degrees) is A </> = 180/M. One can then do least-squares (y 2 ) fitting to estimate both a, </>, 
and their la uncertainties (in practice we fit the modulation curves with the Stokes form of 
the distribution, equation 2). We call this a measurement of polarization. In effect, what is 
measured is the modulation amplitude, a, and it is the knowledge of the detector system, 
expressed in terms of /j, that enables this to be converted to a source polarization amplitude 
via the expression a p = a/ /x. 

This procedure is quite amenable to simulation with Monte Carlo techniques, and here 
we present results of such simulations, both in the case with a 0 = 0 and a 0 > 0. For now we 
ignore background considerations and work only in terms of modulation amplitudes, that 
is, the following results can be considered applicable to any detector system, regardless of 
the /i value that characterizes its polarization sensitivity. 

The simulations proceed as follows, for a given set of true parameter values, ao and </> o, 
we compute a number, M mrrll of simulated data sets each of which has N photons. For each 
simulated data set we bin the resulting angles to form a modulation curve and for each the 



Stokes form of the distribution (equation 2) is fitted to determine best-fit values for Q, U, 
and /, and their la uncertainties, ctq, a U: and 07 , respectively. We can use the fitted Stokes 
parameters to express the results in terms of a and <p using the expressions defined above. 
The la uncertainties, a a and a can also be determined by standard error propagation 
methods, which yields, 

= a (( a 7 //) 2 + ( Oq /{ Q 2 + U 2 ))) 1/2 , ( 6 ) 

and 

a 0 (deg) = (180 /tt) (s/lx 2 ) / (2(1 + a: 2 ))) ((oq/Q) 2 + ( a v /U ) 2 ) 1/2 , (7) 

where, x = Q/U , and 07 , ctq, and cry are the standard la uncertainties on the fitted 
quantities, /, Q and U. The procedure can then be repeated with different values of N. 
All the Monte Carlo simulations described here were done using IDL, and uniform deviates 
were obtained with IDL ’s random number generator randomu. We use a least-squares 
fitting routine within IDL developed by Craig Markwardt that is based on MINPACK-1 
(see http://cow.physics.wisc.edu/craigm/idl/fitting.html). In all of the least squares fits all 
parameters are allowed to vary. 


4.1. Results with ao = 0: Minimum Detectable Amplitude (MDA) 

In the case of ao = 0 the distribution of angles, <p , is uniform, which can be easily 
simulated with a random number generator that produces uniform deviates. We obtain 
from this procedure a distribution of best-fit amplitudes, a, from which we can estimate 
the value cq% that has a 1% chance probability of being measured. This is just the familiar 
MDA described above. Figure 2 shows a comparison of the results from the simulations 
(blue square symbols) versus the analytic expression given above, 4.29/y/N (solid line). 
We computed M sim = 10, 000 simulations for these results, and we used 16 position angle 
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bins for the modulation curves. Figure 2 shows that the ao = 0 simulations are in good 
agreement with the analytic result, giving us confidence that our simulation procedures are 
correct. 


4.2. Results with a 0 > 0: Detection of Polarization 

We can now explore a number of issues with regard to detection of polarization. For 
example, how many counts are needed to measure a modulation amplitude to a particular 
precision, and for a given precision in the amplitude measurement, how accurately can the 
position angle be measured? 

To address these questions we perform additional simulations using true distributions 
with specified amplitudes and position angles. For the illustrative examples below we used 
two different amplitudes, both with 0 O — 0, however, we have explored many different 
cases and all the results summarized here are independent of the particular values of ao or 
0o used. The example distributions described below are; S'(0) = 10 + lcos 2 (0 — 0) and 
S'(0) = 10 + 2 cos 2 (0 — 0). These examples correspond to intrinsic modulation amplitudes, 
ao, of 1/21 and 2/22, respectively. For a specific detector these would correspond to 
polarization amplitudes, a p = l/(21/i), and 2/(22 /i). 

To determine the number of counts, N, needed to reach a polarization sensitivity given 
by the MDP value we require that MDP = 4.29/ /i\fN = a p = a/ /a, and thus a = 4.29 /y/N . 
Note that here /a cancels out, and to achieve MDP values at these amplitudes requires 
(21 2 ) * 4.29 2 = 8116.21, and (22 /2) 2 * (4.29) 2 = 2226.90 counts, respectively. Note that this 
may seem a trivial point, but it’s important to keep the terminology as precise as possible. 
Below we will compare the counts needed to reach a certain MDP, N m d p , and the counts 
needed to measure the same amplitude of polarization to a certain precision, N. While N 
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and N mdp will individually depend on /i, the ratio, N/N mdp , cannot depend on /i, that is, it 
is independent of the detector system employed. 

When Oo > 0 we sample random angles from the true distributions by the so-called 
transformation method. We first compute the cumulative distribution of S((f>), then draw 
uniform deviates, x, and find the corresponding value of (p{x) in the cumulative distribution. 
We use the root finder ZBRENT implemented in IDL to solve for <p given x. We can 
thus draw a specific number of events, N, from the true distribution. For each value of 
N we compute a large number M stm of realizations, to each of which is fitted the Stokes 
distribution to determine the best fit values of /, Q, and U (and thus a and </>), and their 
la uncertainties. We can thus simulate the distributions of a and (p for different N. 


4.3. Single Parameter Confidence Regions 

Here we present the results of our simulations in several ways. First, we used the results 
from M sim = 1, 000 realizations for different values of N to find the mean values of a, and <p 
and their 68% confidence ranges. These are derived for each parameter independently of the 
other, that is, they are 1-dimensional (single parameter) ranges. For example, if one were 
interested in asking the question, is a source (or population of sources) polarized, without 
regard to the particular position angle of the electric vector, then the 1-d distribution would 
be appropriate. We explore the joint, 2-dimensional distributions below in §4.4. 

Figure 3 shows an example simulated distribution of measured amplitudes, a, computed 
with N = 20, 000 events and a o = 2/22 (with <p 0 = 0). The left panel shows the differential 
(binned) distribution, and the right panel shows the same results expressed as a cumulative 
distribution. We also computed similar distributions for the measured position angle, (p. 


Procedurally we use the estimated cumulative distributions to identify the mean of 
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each distribution (ie. , for both a and (f>), and the two values a max and a min that enclose 
68% of the distribution. This then gives the la uncertainty, a ald = ( a max — a min )/2 (we 
also compute the corresponding values for the cj) distribution). We can then compute the 
quantity /3m = cWan/cD, id, which can be thought of as the “number of sigmas” of the 
measurement. We now explore the behavior of several quantities as a function of (3m- 

The first is the ratio N/N mdp (see Figure 4), where N is simply the number of events 
(photons) simulated (ie. the number of observed counts in the modulation curve), and 
N mdp = 4.29 2 /oq is just the number of counts that would be required to reach an MDA 
equal to the true amplitude, cto- This ratio can be thought of as the additional observing 
time required to measure the true amplitude to a given significance compared to the time 
needed to reach an MDA equivalent to the true amplitude. The black diamond symbols in 
Figure 4 show the results of simulations using the example distributions described above 
with M S i m = 1, 000 for different values of N. 

One can also present the results in a slightly different but complementary way. In 
the above simulations we have essentially carried out many simulated observations and 
computed the distribution of “observed” values of the amplitude and position angle, but in 
reality an observer may not have the luxury of making such a large number of independent 
observations of a particular source. All any observer can do is to observe some number, 
N, of photons from the source, construct a modulation curve and fit it as we have done in 
the simulations described above. From this procedure we obtain four quantities, a, a a , <% 
and 00. This constitutes a measurement of the polarization parameters, or more simply, a 
measurement of polarization. Here, the uncertainties, a a , and a^, are obtained for a specified 
confidence level and number of degrees of freedom. For example, for the la (68.3%), single 
parameter confidence ranges we would find the change in each parameter that produces 
a Ay 2 = 1 (while allowing the other parameter to vary in the fitting procedure). For a 
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2-dimensional, joint confidence region at the same level of confidence (68.3%) we would find 
the Ay 2 = 2.3 contour in the a - 0 plane. 

For the diamond symbols in Figure 4 above we used the mean values and la 
uncertainties derived from the distributions computed from many simulated observations 
to obtain /3i d , however, one can also use the “measured” quantities from each simulated 
observation to compute /3 0 b s ,id = a/& a ,id- One can then plot the observed quantities for 
each simulated observation, where now N mdp (a ) is computed using the best-fit value for 
a and the formula a = 4.29 /N^ p . We have done this and show the results in Figure 4 
with the colored symbols. That is, the colored points are these “measured” values from 
individual simulations. The red symbols were computed with N = 10,000, a 0 = 2/22, 

4 > 0 = 0 (deg), the blue used N = 20,000, a 0 = 2/22, 0 O = 0 (deg), and the green with 
N = 24,000, «o = 3/25, and 0o = 22.5 (deg). We see that the distribution of “measured” 
points falls along the same relations as that deduced from the simulated distributions of 
a, and % b, as indeed they should since they are sampling the same distributions, this way 
of presenting the results of the simulations makes a more direct connection with actual 
polarimetric observations, as we only plot quantities that one would obtain directly from a 
single observation. To the extent that actual observations are dominated only by poisson 
counting noise and for which the background is small, then they must fall along the relations 
followed by the simulated observations shown in Figure 4. Indeed, the locus of points traced 
out by the “measured” values from individual simulations can be easily approximated by 
simply plotting curves that intersect the entire swarm of simulated points. Doing this we 
find that N/N mdp = (3° bs ld /9.18 (dashed curve in Figure 4) to very good accuracy. 

We also explored simulations where the true amplitude a 0 was set equal to the MDA. 
Figure 5 shows the cumulative distribution of the “measured” values of f3 from several such 
simulations. One can see from Figure 5 that roughly 60% of the time one would “measure” 



13 


the amplitude, a, at the 3 a level, or better. We emphasize that this relation is based on the 
1-d (single parameter) conhdence range for a, and would be appropriate only for the case of 
addressing the question of the detection of a significant modulation amplitude independent 
of the position angle. We now explore the joint, 2-dimensional conhdence regions. 


4.4. Joint Confidence Regions for a and <f> 

For a 2-dimensional conhdence region we need to find the contour in the a - (p plane 
that encloses a specihed fraction of the best-ht pairs. For a la region (68.3% conhdence) 
this is the contour that satishes A y 2 (a, 0) = 2.3, where Ay 2 = y 2 (a, 0) — Xmin( a best, (ftbest)- 
We again use the Stokes decomposition and compute Ay 2 (J, Q , U) on grids of /, Q and U 
around the best-ht values, I best , Qbest, and Ub es t • We can then convert the grids of Stokes 
parameters into the appropriate values of a and 0 and hnd the boundaries of the region 
that satishes Ay 2 < 2.3. We next hnd the maximum and minimum value on the boundary 
for each parameter. We can then define a aM = ( a mox - a min )/2 and a $ = ( 0 max - <t> m in)/ 2, 
where a max , a min , 0 max and 0 min are the maximum and minimum values on the contour of a 
and 0, respectively. Figure 6 shows a pair of la conhdence regions computed in this fashion. 
Results from two simulations are shown, one with a lower amplitude, aoj ow = 2/24, and 
one with a higher amplitude, — 3/25. Both simulations were performed with </> 0 = 25 
deg (these true parameters are marked by the red diamond symbols), and N = 10000. In 
each case a modulation curve was randomly sampled using the true amplitude and position 
angle, and the simulated data were then htted to determine the best-ht values of the 
amplitude and position angle. These points are shown by the green square symbols. The 
shaded areas show the regions of a and (j) around each best-ht pair that satisfy Ay 2 < 2.3. 
The horizontal and vertical dotted lines mark the maximum and minimum values on the 
regions for each parameter. One can see that the size of the conhdence region grows for 
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smaller intrinsic modulation amplitudes (as one would expect), and that the conhdence 
regions are in general not circular. 

The probability distribution described earlier (§3) can also be used to derive an 
analytic expression for the 2-d conhdence contour in the a - 0 plane for any desired level of 
conhdence. This approach has been investigated by Weisskopf et al. (2010). They derive a 
pair of parametric relations for the values of a and 0 on any conhdence contour (see their 
equation 8). These expressions predict the correct range (extremes) in the amplitude, a, 
but overpredict the range in 0 by a factor of 2, apparently because their original derivation 
neglected a switch from phase angle to position angle (Weisskopf, private communication). 
Thus, if one replaces all occurrences of 0, 0 O and 0 with 2x the respective angle (eg., 
sin 0o > sin 20o ) beginning at their equation (6), one obtains the following expressions; 

a = (oq + A a? c + 2a 0 Aac cos 2(0 — 0 O )) (8) 

and 

0 = 1/2 tan -1 ((a 0 sin 20o + Aac sin 20)/(a o cos 20 o + Aaccos20)) , (9) 

where A ac = (— 41n(l — C)/N) 1 / 2 , with C and N being the desired conhdence level and 
number of detected photons, respectively, and 0 is just a parametric angle that varies 
around the contour. The thick blue curves in Figure 5 were drawn using these expressions 
with N = 10,000, C = 0.683 (ie., ler), and the pair (ao,0o) given by the appropriate best 
ht values (green square symbols). These contour curves provide an excellent match to the 
boundaries of the shaded regions, indicating that the Monte Carlo simulations and analytic 
calculations are in excellent agreement. 

We can now compute the value 0 2 d = abest/&a, 2 d for each particular simulation. This 
again quantihes the “number of sigmas” of the measurement, but now rehects the fact 
that it is a joint measurement of both a and 0 together. Results from a number of such 
simulations are also shown in Figure 4, where we plot the same figure of merit, N/N mc ip, as 
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before, but now using f3 = cibest/ a a ; 2 d- This curve is again quadratic in f3 but rises more 
steeply than the Id relation, because a a2 d is larger than a^id- In agreement with Weisskopf 
et al. (2010) we find that the 2d relation is very well approximated as N/N m d p = /3 2 /4.1, 
and this is the dashed line running through the square symbols. 

We also show in Figure 7 the position angle uncertainty, a^ in degrees (at la confidence) 
as a function of f3 = ab es t/^a,2d- We find that this relation can be very well approximated 
as a^ = 28.5(deg)//3. Thus, a joint 3cr measurement constrains the position angle to better 
than 10 degrees. 


5. Discussion and Summary 

The MDP is a very commonly employed figure of merit to describe the polarization 
sensitivity of a detector system. In their recent paper Weisskopf et al. (2010) have argued 
that “more counts would be needed” to measure the polarization corresponding to the 99% 
confidence MDA rather than to just establish the same level. As we have shown this is 
correct in the sense of a “polarization measurement” being a joint measurement of both 
the amplitude a and position angle <fi. This seems sensible, since it adds an additional 
requirement, that the measured a and (j) both fall within a 2d confidence region. We have 
also confirmed that the additional number of counts required is given by a factor of ss 2.2. 

However, if one were to ask a different question, and were interested in establishing 
simply that a source is polarized, with say, a 3 a measurement of the amplitude, a, without 
concern for the position angle (j), then the Id (one interesting parameter) case is appropriate, 
and no extra counts are needed to measure the MDA. So, in some sense the answer one 
obtains depends on the question asked. 

The MDP value has often been used to estimate observing times required to reach 



particular sensitivity levels. The results shown here demonstrate that exposure requests 
should clearly match the measurement goals of the desired scientific program. For example, 
if source modeling requires a joint measurement of both polarization parameters, then the 
appropriate time to reach a required precision for two interesting parameters should be 
requested. 
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Fig. 1. — Geometry relevant for photoelectric polarimeter measurements. Photons travel 
from figure top to bottom. Photoelectrons are preferentially emitted in the plane perpendic- 
ular to the photon direction of travel (in this case the X - Y plane), and their initial direction 
is measured by the azimuthal angle 0. 
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Fig. 2. — The amplitude, a\% = MDP(n = 1), vs the number of counts, N , obtained from 
simulations described in §4.1 (blue squares), and the analytic formula (solid curve, §3). 
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Fig. 3. — Example Id distributions in the amplitude, a, computed from Monte Carlo simula- 
tions. The left panel shows the differential distribution (binned), and the right panel shows 
the cumulative distribution. To estimate the mean we find the midpoint (triangle sym- 
bol), and the la extremes, a m i /n and a ma x at 0.165 and 0.835, respectively (square symbols, 
enclosing 68% of the distribution). 
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Fig. 4. — Plot of N/N mdp (a ) as a function of /3, the “number of sigmas” of the measurement. 
Both the Id amplitude distribution (cr a ld , independent of the position angle) and the 2d 
joint distribution (cr ai 2 d), are shown. The black diamond symbols are derived from the 
results of M sim = 1000 simulations for different values of N (see discussion in §4.3). The 
colored symbols are the results of individual simulations where /3 is derived from the best- 
fit amplitude, a, and its Id, la uncertainty, cr a ,id- The red symbols were computed with 
N = 10,000, ao = 2/22, 0 o = 0 (deg), the blue used N = 20,000, ao = 2/22, 0 O = 0 (deg), 
and the green with N = 24,000, a 0 = 3/25, and 0 O = 22.5 (deg). The solid dashed curve 
for the Id case is given by N/N m d p (a ) = /3 2 / 9.2. The black square symbols show the 2d 
confidence region results (see discussion in §4.4). The curve running through the 2d results 
is given by N/N mdp (a ) = /3 2 /4.1 
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Fig. 5. — Plot of the cumulative distribution of j3 = a/cr aj id for several different simulations 
all satisfying the condition that a = MDA. A vertical line is plotted at /3 = 3 (the nominal 
3cr detection criterion), and which is close to the most probable value (the distribution is 
not exactly gaussian, ie. symmetric). See the discussion in §4.3 for more details. 
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Fig. 6. — Confidence regions and contours in the (a, 0) plane. Results from two simulations 
are shown. Simulated modulation curves were computed with N=10,000 counts for two 
different intrinsic amplitudes both with a position angle of 25 degrees. These values are 
marked by the red diamond symbols. The resulting best fit parameter values (green squares) 
and confidence regions (shaded areas) are shown. The shaded regions are the Ay 2 < 2.3 
regions (lcr for 2d confidence regions). The horizontal and vertical dotted lines denote the 
extremes in each parameter and are used to compute a at 2 d and cr^d- The blue contour curves 
were computed from the analytical expressions in §4.4. 




